Log(l/x) Gluon Distribution and Structure Functions 
in the Loop-Loop Correlation Model 
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We consider the interaction of the partonic fluctuation of a scalar "photon" with an external color 
field to calculate the leading and next-to-leading order gluon distribution of the proton following the 
work done by Dosch-Hebecker-Metz- Pirner [lj. We relate these gluon distributions to the short and 
long distance behavior of the cross section of an adjoint dipole scattering off a proton. The leading 
order result is a constant while the next-to-leading order result shows a ln(l/a;) enhancement at 
small x. To get numerical results for the gluon distributions at the initial scale Q% = 1.8 GeV 2 , 
we compute the adjoint dipole-proton cross section in the loop-loop correlation model. Quark 
distributions at the same initial scale are parametrized according to Regge theory. We evolve quark 
and gluon distributions to higher Q 2 values using the DGLAP equation and compute charm and 
proton structure functions in the small- a; region for different Q 2 values. 
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I. INTRODUCTION 

The understanding of deep inelastic scattering (DIS) in 
the small- a; regime remains one of the challenges in quan- 
tum chromodynamics (QCD). Perturbative and nonper- 
turbative physics are important for a complete picture of 
the small- a; limit of structure functions. In this work, we 
combine perturbative and nonperturbative approaches to 
describe charm and proton structure functions (or quark 
and gluon distributions) in the small- a; region. 

The basic idea is the description of high-energy scatter- 
ing in QCD by studying the eikonalized interaction of en- 
ergetic partons with external color fields. In ref. 0, this 
semiclassical method has been compared with the parton 
model to extract the leading order (LO) and next-to- 
leading order (NLO) gluon distribution. The main idea 
is simple and goes back to Mueller 



One calculates 

the gluon production cross section for a scalar "photon" , 
which directly couples to the gluon field, in an external 
color field. The calculation of one-gluon production leads 
to the LO gluon distribution and the one of two-gluon 
production to the NLO gluon distribution. The LO result 
turns out to be a constant, xg^(x,Q 2 ) oc const., char- 
acterising the averaged local field strength in the proton. 
The NLO result shows a logarithmic increase at small 
x, xg^(x,Q 2 ) cx ln(l/a;), and is sensitive to the large 
distance structure of the proton. 
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In order to make numerical estimates for the gluon dis- 
tributions, we relate the LO and NLO gluon distributions 
to the scattering of a dipole in the adjoint representation 
of 5/7(3) on a proton. Then, we use the loop-loop corre- 
lation model (LLCM) |3j to compute this adjoint dipole - 
proton cross section at the center of mass (cm.) energy 
of y/so « 20 GeV. The two gluons emerging from the 
scalar "photon" represent the adjoint dipole and a fun- 
damental quark-diquark dipole models the proton in the 
LLCM. The correlation between the two dipoles given 
by Wegner- Wilson loops is evaluated in the Gaussian 
approximation of gluon field strengths. The perturba- 
tive interactions are obtained from perturbation theory 
and the nonperturbative gluon field strength correlator 
is parametrized in line with simulations o in the lattice 
QCD. The calculation of scattering cross sections in the 
loop-loop correlation model has been quite successful at 
low energies ^ « 20 GeV 0, 0, S B . 

For small virtuality Q 2 , the perturbative wave function 
of the two gluons emerging from the scalar "photon" is 
unrealistically extended at the cndpoints of the longitudi- 
nal momentum fraction a = 0, 1. This problem is similar 
to the quark-antiquark wave function of the transverse 
photon. The phenomenological solution we propose is to 
give the gluon a constituent mass rog of the order of the 
rho mass. This mass modifies the wave function at small 
transverse gluon momenta k± ensuring "confinement" for 
the two gluons, but does not affect the perturbative part 
of the two-gluon wave function at high k±. 

We calculate the x-dependence of the gluon distribu- 
tion at a scale Q 2 , = 1.8 GeV 2 which corresponds to 
the upper limit where the nonperturbative input is still 
credible. The quark distributions at the same scale are 



parametrized in line with Regge theory |8(. We evolve 
both distributions to higher values of Q 2 using DGLAP 
equations Q. With quark and gluon distributions, we 
compute charm and proton structure functions at small 
x for different Q 2 values in good agreement with experi- 
mental data. 

The outline of the paper is as follows: In section [Q] 
we review the formulas for the calculation of the LO and 
NLO gluon distribution following Q], introduce a gluon 
mass into the formalism and rewrite the final expressions 
for gluon distributions in terms of adjoint dipole-proton 
cross sections. We compute the adjoint dipole-proton 
cross section in the loop-loop correlation model in sec- 
tion HTU In section |^ we determine the x dependence of 
the quark distributions at the initial scale Qq according 
to Regge theory. The intial gluon and quark distributions 
are evolved to higher virtualities Q 2 using the DGLAP 
equation. In section we present the results for the 
charm and proton structure function versus Bjorken-ir. at 
different Q 2 values. Finally, in section IVTl we summarize 
our results. 
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FIG. 1: The process \ ~~ > 9 in an external color field: 
(a) one gluon exchange, (b) resummed gluon exchange. 



x + = xq + X3 the light-cone variable, d± = d/db±, 
A A = A b t b the external color field, Ca = N c the Casimir 
operator in the adjoint representation, N c the number 
of colors and the polarisation vector of the outgoing 
gluon. 

Resumming the gluon exchange to all orders as shown 
in Fig. ^b) , the gluon initially created at the vertex \99 
aquires a non-abelian eikonal factor on the way through 
the external color field 



II. THE SEMI-CLASSICAL GLUON 
DISTRIBUTION 



U, 



(b±)=P exp 



dx + A A (x+,b x ) 



In this section we review the formulas for the compu- 
tation of the LO and NLO gluon distribution [3|| in the 
semiclassical approach and parton model following the 
work by Dosch, Hebecker, Metz and Pirner [lj. In addi- 
tion, we introduce a gluon mass in the formalism to take 
into account the spatial localization of the two gluons 
due to confinement. We derive the modifications due to 
massive gluons and rewrite the final results for the gluon 
distributions in terms of adjoint dipole-proton cross sec- 
tions. 



The path ordering along the way x+ is denoted by P. 
With 

W b A Jr A ) = U A (b ± - r A /2)U A \b x + r A /2) - 11 , (1) 

where the eikonal factors U A (b± — r A /2) = U A _ QO (b± — 
r A /2) and U A ^ (b±+r A /2) come from the scattering of two 
gluons (moving in opposite directions) with transverse 
distance r A off the proton, the semiclassical (sc) cross- 
section for gluon production at leading order becomes 
(cf. ref. 0) 



A. Gluon distribution at leading order 

In order to extract the gluon distribution, we consider 
a scalar field \ ( or scalar "photon") coupled directly to 
the gluon field through the interaction lagrangian 



X 



rX tr (F^F 



fJ,U\ 



with F M „ = F^t a , the gluon field strength F* v = 
[d^A a v - d u A° ~ gf abc A h ^Al\ , the gluon field A a and the 
SU(N C ) group generators in the adjoint representation 
t a . In the high-energy limit, the scattering of a x-particle 
off an external color field A a as shown in Fig. [Tfa), has 
the following amplitude 

T a (b±) =^ J dx +^ [t a (e* ± d ± )A A (x + ,b ± )] ■ 

Here b± denotes the impact parameter in transverse 
space, q the four momentum vector of the x - P ar ti c l e i 



r(°) = 



A 2 



d 2 b { 



[dr A W A (r A )} 



r A =0 



(2) 



In the following an average over the external color fields 
underlying the quantity W A ^ is implicitly understood. 

On the other hand, the leading-order parton model 
(pm) cross section for the partonic process \9 ~ * 9 shown 
in Fig. |21 with the proton described by the gluon distri- 



FIG. 2: The process X9 — * 9 in the parton model. 



bution is given by 



r(°) 



xgM(x,Q 2 ). 



(3) 



Identifying the semiclassical cross-section with the 
partonic one one finds the leading-order gluon dis- 
tribution 



xg<-°\x,Q 2 ) 



1 



1 



2ir 2 a s 2C A 



d 2 b t 



[d rA W b A ± (r A )] 



(4) 

Using the expression for the adjoint dipole (d)_^ proton 

(5) 



) - HJ 1 

(p) cross section in the semiclassical approach P, [l(| 



°%<Ta) = -\ J d 2 b±_ txW b A ± (r A ) 
and the identity 



d 2 b± 



[dr A W A {T A )] 



r A =0 



2 d 2 b ± tvW b A (r A ) 



one can rewrite eq. (@J in the useful form 



xg(°\x,Q 2 ) 



4ir 2 a s 2C A 



%-[dlAl{r A )] r „ 



r A =0 

(6) 



(7) 



This equation shows that the LO gluon distribution de- 
pends on the small distance behavior of the adjoint 
dipole-proton cross section. To get numerical results for 
the gluon distribution, we use in the next section a model 
to compute the adjoint dipole-proton cross section. 



B. Gluon distribution at next-to-leading order 





FIG. 4: Relevant contributions in the NLO semiclassical cal- 
culation. 



where a and I — a are the longitudinal momentum carried 
by the two gluons, k± and k'j_ are the transverse momenta 
of one of the two gluons before and after the interaction 
with the external field, TV 2 = a(l — a)Q 2 and W^(k±) 
is the Fourier transform of W^[rX) given in eq. Q. 

In equation JSJ we have introduced a gluon mass tuq 
into the gluon propagator. This mass mimics the con- 
finement of the two gluons emerging from the x _ P ar ticle 
and modifies their perturbative wave function in the non- 
perturbative, small k± region for small virtualities Q 2 or 
momentum fractions a = 0,1. Phenomcnologically, the 
gluon mass has a similar effect as a constituent quark 
mass in the perturbative quark-antiquark wave function 
of the transverse photon 11]. We use for the gluon mass 
mo = 770 MeV, i.e., approximately half of the glueball 
mass. 

To compute the NLO contribution to the gluon dis- 
tribution in the parton model, we have to calculate the 
diagrams shown in Fig. [5] With z — Q 2 /(Q 2 + s) = x/y, 



For the semiclassical calculation of the gluon distribu- 
tion at next-to-leading order, the three diagrams shown 
in Fig. |3 are relevant in the high energy limit. Resum- 
ming the interaction with the external field to all or- 
ders, i. e., repeating the step leading from Fig. Ufa) to 
Fig-H^b), we obtain the result shown in Fig. In the 
first diagram the incoming ^-particle splits into two glu- 
ons before interacting with the target. The subsequent 
scattering of the two gluons off the exteral color field 
is treated in the eikonal approximation. In the second 
diagram the two fast gluons are created through a xggg 
vertex in the space-time region of the external color field. 




FIG. 5: Diagrams contributing to the parton model at NLO. 

where y is the fraction of the target momentum carried 
by the struck gluon, the parton model cross-section at 
NLO reads 




FIG. 3: Diagrams contributing to the \ — > gg amplitude in 
the high-energy limit. 

One can show pj that the amplitudes in Fig.0|produce 
the following NLO contribution to the cross-section 



A 2 



da 



32(2tt) 6 J a(l-a) 



dk'; 



d 2 b. 



(8) 



d 2 k t 



N 2 Sij + 2kikj 
N 2 + k 2 



W A (k' ± - fc x ) 



G 



<r P m(x, Q 2 ) = [ - U 0) <5(1 -x) + <7«(x, Q 2 )} yg b (y), 

J x Z L J 

(9) 

with = (7rA^)/4 and the d = (4 + s) dimensional 
coupling \d = A/i~ e / 2 . Since we are interested in the 
high energy limit, we extract only the leading term in 
ln(l/x) from the expression in ©. We regularize the 
bare gluon distribution g b (x) in the MS scheme 



9b(x) = g{x^i 2 ) 



- + IE — ln(4-7r) 



g(g,v 2 ), 



with P gg denoting the gluon-gluon splitting function. 
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The parton model cross-section at NLO finally becomes 



_ (o) 



1 dz 



(10) 



6(1 -z) + 



2tt 



where 



Cf s (z) = P 53 (z)ln 



1 - z 



UC A 



z J 6z(l — z) 



Expanding the full gluon distribution up to the order 
ln(l/x) 

xg(x, y 2 ) = xg^ (x, fi 2 ) + xg {1) (x, ^ 2 ), 
and identifying ijSJ with (fTU|) . we obtain 



da 



d 2 k\ 



4(2tt) 7 J a(l - a) 
N 8jj + 2kjkj T ~ rA , 



dk'l / d 2 b± 



N 2 + k\ + m% 



W b A ± (k' ± -k x ) 



P 9g (z)ln(^)+C g MS (,) 



2/5 ( °W), 
(11) 



where g^\x,[i 2 ) is given by 10}. When evaluating equa- 
tion we keep only ln(l/x) terms, i.e., we use 



Cf^(z) 



2Ca 

z 

z 



In 



1\ 11 



z 12 



Using the variable z = Q 2 /(Q 2 + M 2 ), with M 2 = 
k l2 /(a(l — a)), and introducing a scale k 2 such that 
Aq C£) <C k 2 <C Q 2 , the first contribution in lfTT|) can 
be split into a hard and a soft component Q, 0] 



da 



4(2tt) 7 J a(l-a) 

AT 2 ^,- + 2fcife 



dfc^ / d 2 6. 



d 2 fc. 



A 2 + fc^ + : 



4tt 3 f x z 1U (zk 2 ) &± 

+- f 1 - r dk> 2 f(k<i), 

n Jx z Jo 



(12) 



where 
f(k'l) 



(2ir) 2 r 2 J (2n) 2 r A 2 



tr 



H = 



{rA ' r '* )2 [a 2 K (a)} \b 2 K (b) 



r A r A 



2(r A • r' A ) 



l ^2 



{[aKi{a)\ [bKAb) 
+ [aKi(a)} \b 2 K Q (b) 
+ [a 2 K (a)] [SiTi(S)]}, 



with a = vhqt a and b = m,Qr' A . Note that the gluon 
mass does not influence the hard part which is calculated 
at leading twist. The lower bound of the z integration, 
x < z, is a kinematical limit ensuring that the invariant 
mass of the two produced gluons cannot be larger than 
the total center-of-mass energy available 0, . 

Inserting (|12J) in (|ll|l . the ln 2 (l/x) terms from the 
semiclassical and the parton calculations cancel, so that 



xg^\x,fi 2 ) = In 



1 

X I | Ti 

a 



dk' 2 f(k'j 



7T 



C A 



In 



11 
12 



xg 



(0) 



For k 2 — fi 2 exp( 11/12), the second term also drops out 
and one obtains 

^(^^-M- / dk'tfik' 2 ). (13) 

7T \XJ Jo 

After the integration over k'l and rV the NLO correction 
to the gluon distribution becomes "0, 0] 



x / d 6j.tr 



^W< T (^) (14) 



where 



F(z) 



4eA^ 27E 

F 2 (z) + zK^z)Ki(z) + \K 2 {z) 
7 ^K 2 2 (z) - zK^z)K 2 {z) + Kl{z). 



In the large- N c limit [T2|, one obtains the relation 



d 2 6_Ltr 



W£(r A )w£(r A ) = -2 I d 2 b ± tr W&(r A ) 
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which allows us to rewrite the NLO correction in terms 
of the dipole-proton cross section JSJ 



:(^ 2 ) 



-^F{m G r A ) *%(r A ). 

(15) 



The NLO gluon distribution depends on the large dis- 
tance behavior of the adjoint dipole-proton cross sec- 
tion. For small rj^ the function F(rriGi'X) reduces to 
l/( m G r ^) 2 ! an d eq. (15) shows qualitative agreement 
with eqs. (49)- (50) of Mueller 0. To obtain numerical 
results for the LO and NLO gluon distribution, we cal- 
culate the dipole-proton cross section at some low cm. 
energy ^/Jq within the loop-loop correlation model in the 
following section. 



III. DIPOLE-PROTON CROSS SECTION FROM 
THE LOOP-LOOP CORRELATION MODEL 

In this section we use the loop-loop correlation model 
(LLCM) to compute the cross section of an adjoint 
dipole scattering off a proton, <r^(r A ,so), at the cm. 
energy ^/so = 20GeV. Then, we insert the resulting 
dipole-proton cross section in eqs. and fT5|l to calcu- 
late the LO and NLO gluon distributions. 

The loop-loop correlation model is based on the func- 
tional integral approach to high-energy scattering in the 
eikonal approximation 0, 0, El • Its central elements 
are gauge-invariant Wegner- Wilson loops. With a phe- 
nomenological energy dependence, we have shown that 
the LLCM allows a unified description of high-energy 
hadron-hadron, photon-hadron, and photon-photon reac- 
tions @,E3>0] and of static properties of hadrons [l5llTf^ | 
in agreement with experimental and lattice data. In this 
work, we do not need a phenomenological parametriza- 
tion for the energy dependence since this is generated by 
the formalism outlined in the sections before. 

In the framework of LLCM the adjoint dipole - proton 
cross section reads 0, 0] 

a LLCM( r A> s o) = 2 J d2b± f~~^~ J ^ r Fdz g \ip p (z q ,j%)\ 2 
x(l-S A ^(s ,b ± ,r Al z qi r^) . (16) 



Here the correlation of two light-like Wegner- Wilson 
loops 



S A ^(s ,b x ,f A ,z q ,f r ) = — ^ '-G 



w A [c g ]) G {wr[c q ] /i; 



(17) 

describes the elastic scattering of a dipole in the funda- 
mental (T) with a dipole in the adjoint (A) representa- 
tion of SU(N C ). In the present case where a x-particle 
scatters off a proton, the adjoint color-dipole is given by 



the two gluons in the color-singlet state emerging from 
the x-particle and the fundamental color-dipole is given 
in a simplified picture by a quark and diquark in the pro- 
ton [3^| . Each color-dipole is represented by a light-like 
Wegner- Wilson loop [l8| 



W r [C] = fr r -Pexp 



-ig<b dz»A a Jz)t a r 



where the subscript r indicates the representation of 
SU(N C ), fr r = Tr r (- • ■ )/Trll r is the normalized trace 
in the corresponding color-space with unit element H r , 
and Afj_(z) — A"(z)£° represents the gluon field with the 
SU(N C ) group generators in the corresponding represen- 
tation, t", that demand the path ordering indicated by 
V on the closed path C in space-time. Physically, the 
Wegner- Wilson loop represents the phase factor acquired 
by a color-charge in the SU(N C ) representation r along 
the light-like trajectory C in the gluon background field. 

The color-dipoles have transverse size and orientation 
r A jr. The longitudinal momentum fraction of the dipole 
carried by the quark (gluon) is Zn (%)■ The impact pa- 
rameter between the dipoles is [ljj 

b± = Tg + (1 Zg)r A Tq (1 Zq)r F = r Acm ?> cm , 

where f qi f qq , r g , and r q are the transverse positions of 
the quark, diquark, gluon and the gluon moving in the 



opposite direction, respectively. With 7> 



r a and 



f g: the center of light-cone momenta of the 



two dipoles are given by r, 



Acm — z g T g + (1 z g) r g an< I 
T^cm = z q f q + (1 — z q )r qq . Figure |S| illustrates the (a) 
space-time and (b) transverse arrangement of the dipoles. 

The QCD vacuum expectation value (. . .)g in the loop- 
loop correlation function (|17|) represents functional inte- 
grals over gluon field configurations: the functional 
integration over the fermion fields has already been car- 
ried out as indicated by the subscript G. The model we 
use for the QCD vacuum describes only gluon dynam- 
ics and, thus, implies the quenched approximation that 
does not allow string breaking through dynamical quark- 
antiquark production. 

The r> and z q distribution of the fundamental color- 
dipole in the proton is given by the proton wave function 
ipp. We use for the proton wave function the phe nomeno- 



logical Gaussian Wirbel-Stech-Bauer ansatz 



prion 

m 



tp p (z qi f^) 



2ttS^N p 



-(z„ 



7(4A 2 



l^| 2 /(4S=) 



The constant N p is fixed by the normalization of the wave 
function to unity, the extension parameter S p is approxi- 
mately given by the electromagnetic radius of the proton 
and the width Az p = w/(\^2m p ) [2f| is determined by 
the proton mass m p and the value w = 0.35 — 0.5 GeV 
extracted from experimental data. We adopt the values 
Az p = 0.3 and S p = 0.86 fm which have allowed a good 
description of many high-energy scattering data in our 
previous work 0. 
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(b) 



C 3 = dS 9 




FIG. 6: High-energy scattering of a fundamental dipole with 
an adjoint dipole in the eikonal approximation represented 
by Wegner- Wilson loops: (a) space-time and (b) transverse 
arrangement of the Wegner-Wilson loops. The shaded areas 
represent the strings extending from the quark (gluon) to the 
diquark (antigluon) path in each color dipole. The thin tube 
allows to compare the field strengths in surface S q with the 
field strengths in surface S g . The impact parameter b± con- 
nects the centers of light-cone momenta of the dipoles. 



field strengths. These steeps lead to the following result 



G 



1 



W*[C B ]) a (wr[C q ] /c 



AT|_i 

N c + 2 
2(N C + 1) 

N c ~2 
2(N C -1) 



exp 



.N e 



exp 



exp 



l 2 X 



r l 

*2* 



(18) 



which, for N c — 3, corresponds to the well-known SU(3) 
decomposition 

3^8 = 3® 15 9 6. 

The function \ '■= X P + Xnc + Xc P has the following 
form 



X P = 9 2 (r g - 


f^iD'P (f g - 




+ 9 2 (r§- 


f qq ) iD' p {2) {r- g - 


- r qq 


~ 9 (r g - 


r qq ) iD' p {2) (f g - 


- T qq 


~ 9 (r- g - 


r q ) l Df ) (f g - 




NP 7r 2 G 2 (l- 
3(iV 2 - 







(19) 



V( 2 ) 



' qq I 



iD'V\r g -r qq )-iD'f\?- g -r q ) 



NP _ K G2K 

Xc ~ 3(N 2 - 1) [Ta ' Tt) 



(20) 
(21) 



x / dv A / dvr iD^ (f g + v A f A - f q - Vrfjr) , 
Jo Jo 

with the perturbative {iD p ) and nonpcrturbative 
(iD 1 ^ and iD^(z±)) correlation functions in transverse 
space 



iD' p {2) (z ± ) = j-K (m G \zl\), 



1 



(22) 



^^' 1 (2) (4) =Tra 4 [3 + 3(\z ± \/a) + (\z ± \/a) 2 ] exp(-|^|/a), 

(23) 



The computation of the loop-loop correlation for one 
dipole in the fundamental and the other one in the adjoint 
representation of SU(Nc) can be found in detail in 0, 
Il5l Ha |. The main steeps of this computation are the 
transformation of the line integrals in the Wegner-Wilson 
loops into surface integrals with the non-Abelian Stokes 
theorem, a matrix cumulant expansion, and the Gaussian 
approximation of the functional integrals in the gluon 



iDW(z ± ) = 2ir a 2 [1 + (\z ± \/a)] exp(-|4|/a) 



(24) 



We have introduced in the perturbative component \ P 
the same effective gluon mass mc = 0.77 GeV as before 
to limit the range of the perturbative interaction in the 
infrared region and a parameter M 2 — 1.04 GeV 2 which 
freezes the running coupling in the quenched approxima- 
tion at the value g 2 '(zi) / '(4tt) = a s (kl = 0) = 0.4 0, 



9 ( z± 



48ir 2 



(33-2N f ) In (|4 



M 2 )/A 2 



QCD 



(25) 
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In Eq. ij^ K denotes the th modified Bessel func- 
tion (McDonald function). The non-perturbative cor- 
relators (|2f)|> and (|21|) involve the gluon condensate 



G 2 := (^^(0)F^(0)) = 0.074 GeV 4 , the param- 
eter k = 0.74 that determines the relative weight of 
the two different components and the correlation length 
a = 0.302 fm that enters through the non-perturbative 
correlation functions D and D%. 

The component x P describes the perturbative interac- 
tion of the quark and diquark of the dipole in the proton 
with the two gluons of the adjoint dipole emerging from 
the x-particle as evident from i|19|) and Fig.HJ). The com- 
ponent Xnc nas the same structure as \ F an d gives the 
non-perturbative interaction between the quarks and glu- 
ons of the two dipoles. The component \c P shows a dif- 
ferent structure: the integrations over v A and v F sum non- 
perturbative interactions between the strings (flux tubes) 
confining the quark and diquark or the two gluons in the 
dipoles as visualized in Fig. HJ>. As shown in [l5l . l2l"l |22| . 
the Xc P component leads to color-confinement due to 
a flux tube formation between a static quark-antiquark 
pair. Manifestations of confinement in high-energy scat- 
tering have been analysed in [l4|. 

Since the squared proton wave function \ip p (z q , r>)| 2 
is invariant and the %-function changes sign under the 
replacement (7> —> — 7>, z q —* 1 — z q ) 

X(h,Zg,f A ,l- Z q ,-J%) 



-x{b±,z g ,r Al z qi rr) 



only the real part of the exponentials in eq. I|18l) survives 
in the integration over f> and z q so that one obtains 



W A [C g ]W^[C q 



1 



JVf- 



1 

-2 



Nc ■ 
2 X 



2(N C 
N r - 



■ cos 



2(N C -1 



■ cos 



2 X 



2 X 



(26) 



The above expression describes multiple gluonic interac- 
tions between two dipoles since (x P ) 2 represents the per- 
turbatively well-known two-gluon exchange and (x^) 2 
the non-perturbative two-point interaction in the dipole- 
dipole scattering 0] . The higher order terms in the ex- 
pansion of the cosine functions ensure the S'-matrix uni- 
tarity condition which becomes important at very high 
cm. energies ffi Il3l l23| . 

The adjoint dipole-proton cross section obtained with 
the above ingredients at the cm. energy ^/so ~ 20 GeV 
shows color-transparency for small dipole sizes, 



ct *cm( s o>»a) fa 9.6 r 2 



(27) 



and linear confining behavior at large dipole sizes, 

o'liOM^o,^) oc \r A \ ■ (28) 



Inserting eqs. i|27|) and l|25[) into J7J, one obtaines for the 
LO gluon distribution Q at virtuality Q\ = 1.8 GeV 2 



xg^(x,Q 2 ) 



4n 2 a s (Q 2 ) 2C A 



2 ■ 9.6 = 0.81 . (29) 



With our result for cr^ CM (s , r A ), the NLO gluon dis- 
tribution i|15|) at the same virtuality reads 

xg (1 \x,Ql) = 0.89 In (^) . (30) 



IV. DGLAP EVOLUTION AT HIGH Q 2 

In this section we give the gluon and quark distribu- 
tions at an initial scale Qq. Their evolution to higher 
values of Q 2 is obtained by the DGLAP equation 



Q 2 dr 



Qi{x,Q 2 ) 



(31) 




Pq iB \ 








. \g^Q 2 ) 



with the splitting functions P xy being at leading order 
Q 2 independent. We use the resulting gluon and quark 
distributions to compute the charm and proton structure 
function at different x and Q 2 values. 

The gluon distribution computed in the previous sec- 
tion reads 



xg{x,Ql) = A[l + B\n(l/x)} 



(32) 



with A = 0.81 and B = 1.1 for Q 2 = 1.8 GeV 2 . Our re- 
sult, however, is only expected to be valid at low x values. 
For large x values, x > xqrv = 0.15, we use the Gluck- 
Reya-Vogt (GRV) gluon distribution 0. To match to 
this gluon distribution at x = sqrv , we introduce a scale 
Xq in our gluon distribution 



G(x,Ql) = xg(x, 



A[l + BJn(x /x)] (33) 



which takes into account the neglected constant term in 
the NLO calculation where only the leading ln(l/x) terms 
have been kept. For Qg = 1.8 GeV 2 and xqrv = 0.15, 
we obtain xq — 0.1454. 

To calculate the proton structure function F%(x, Q 2 ), 

F${x)=X J2 e 2 q [q(x)+q(x)], 



with the DGLAP evolution (|31|) . we need quark and 
gluon distributions, since they are coupled to each other. 
For the computation of F%(x, Q 2 ), however, only two lin- 
ear combinations of quark distributions are required 

T = x(u + + c+ + t + ) - x(d + + s + + b + ), 
E = x(u+ + c+ +t+) +x(d+ +s+ +b+), 
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since 



10~ 4 10~ 3 10~ 2 



5E + 3T 
L8 



with q + = q + q and q = u,d, s, c, t, b. Performing linear 
combinations in we can directly check that T, E and 
G evolve according to the following DGLAP equations 



Q 2 d Q 2T(x,Q 2 ) 



2n 



Q 2 dr 



E(x,Q 2 ) 
Q2 \G(x,Q 2 ) 



2tt 



1 x<% (P qq 2n f P q; 



Pa, 



Pa, 



fmQ 2 ) 



which means that T evolves alone, while E is coupled 
with G. 

We determine the x-dependence of the T and S dis- 
tribution at the scale Qq following the ideas from |25j . 
First, the T, E and G distributions are required to have 
a common singularity structure in the complex j-plane 
according to Regge theory. Since the initial gluon distri- 
bution (|33|l has a double pole in j = 1 , consequently, also 
the T and S distribution are required to have a double- 
pole pomeron term. Secondly, we add a reggeon contribu- 
tion (coming from the exchange of meson trajectories ao 
and /) to the quark distribution. In the gluon distribu- 
tion we negelect this term since the reggeon is expected to 
be constituted of quarks. Moreover, we expect that the 
pomeron, having vacuum quantum numbers, does not 
distinguish between quark flavours, i. e., the pomeron 
decouples from the T distribution. In contrast, the E 
distribution contains a pomeron and a reggeon compo- 
nent. The initial distributions, thus, read 



T(x,Ql) = Tx aa {\-xY, 

E(x,Ql) = [a\n(l/x) + b + dx ao ]{\-x) K 



(34) 



with the reggeon intercept ao = 0.4 and the powers a — 
3 and k = 2 of (1 — x) taking into account daughter 
trajectories in Regge theory. For x > xgrv, we again rely 
on the GRV distributions since our T and S distributions 
are only valid at small x. The parameters r and b are 
fixed to ensure continuity between our distributions and 
GRV's ones at x = xqrv- Finally, we are left with 2 
parameters, a and d, which we determine by fitting the 
F2(x,Q 2 ) experimental data. 



RESULTS 



With the previous T, E and G distributions at Qq — 
1.8 GeV 2 as an initial condition in the DGLAP equation, 
we have fitted the F% experimental data within the do- 




FIG. 7: The charm structure function, F2 (i, Q 2 ), as a func- 
tion of the Bjorken- variable x at different virtualities Q 2 . 



Q 2 > Ql = 1.8 GeV 2 

X < £GRV = 0.15 , 



(35) 



cos(0 t ) 



IQ 1 > 49 GeV 2 



2 x m p 2m 2 



Here m p denotes the proton mass, 9 t is the scattering 
angle in the t-channel and cos(# t ) has been extended to 
the whole complex plane. The last condition of cos(# t ) 
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FIG. 8: The proton structure function, FJ^a;, Q 2 ), as a function of the Bj or ken- variable x at low virtualities Q 2 



taken from |26J ensures that Regge theory is applica- 
ble. We have used 922 experimental points coming from 
HI mUlJIII, ZEUS HlUni BCDMS p3, E665 jH 
and NMC |34| to adjust our parameters within the given 
region. We have obtained the following values 

a = 0.45933 ±0.00425 , 

b = —1.9598 (fixed by continuity) , 

d = 6.0408 ±0.0456 , 

t = 0.55628 (fixed by continuity) . 



In Figs. [SJ IH1 and ^| we show the results for the 
proton structure function, F^ix, Q 2 ), as a function of 
the Bjorken- variable x between the virtualities Q 2 = 



1.8 GeV 2 and Q 2 = 3000 GeV 2 . Our results are in good 
agreement with experimental data for all x and Q 2 val- 
ues. 



To test our gluon distribution, we compute the charm 
structure function which depends on the gluon distribu- 
tion as follows |33 



F 2 c (x,Q 2 ) = 2ei 



a s (Q 2 + 4m 2 c ) 
2tt 

d£g(tQ 2 + 4m 2 c )f(x/(i,Q 2 ) , 
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FIG. 10: The proton structure function, F2(x,Q 2 ), as a function of the Bjorken-variable x at high virtualities Q 2 . 
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with 



[A-^)x\l-x)-- 



'X 

.2 
2m? 



x\\ - x) + plx' 2 (1 - Zx) - [1 2 Z 3 



= \ l- 



2xfi 



L 



l-x ' 
1 + w x 



log , ., 

1 — V 

l + 2fi 



We have adopted a value of 1.25 GeV for the charm 
quark mass. The predictions for F%(x, Q 2 ) obtained from 
our model arepresented together with the experimental 
HERA data [Ifl S H3 in Fig. □ The good agreement 
of our predictions for F^iXjQ 2 ) with experimental data 
clearly shows that we obtain a resonable gluon distribu- 
tion. 



VI. CONCLUSIONS 

We have considered the interaction of a scalar "pho- 
ton" , which directly couples to the gluons, with an exter- 
nal color field to extract the leading and next-to-leading 
order gluon distribution. We have closely followed the 
previous work by Dosch, Hebecker, Metz and Pirner 
in which the semiclassical approach (where the partonic 
fluctuations of the "photon" interact with the proton in 
the eikonal approximation) has been compared with the 



parton model to get the gluon distribution of the proton. 
The leading order result is a constant while the next-to- 
leading order result has a ln(l/x) rise at small x. 

We have been able to relate the leading and next-to- 
leading order gluon distribution to the short and long 
distance behavior of the cross section of a dipole in the 
adjoint representation of SU(3) scattering off a proton, 
respectively. In addition, a gluon mass has been intro- 
duced to take into account nonperturbative effects in the 
small-fci region of the perturbatively derived wave func- 
tion of the two gluons emerging from the scalar "photon" , 
in analogy to the constituent quark mass in the quark- 
antiquark wave function. We have computed the adjoint 
dipole-proton cross section in the loop-loop correlation 
model [3j to obtain numerical results for the gluon dis- 
tributions at the initial scale Qq — 1.8 GeV 2 . 

Quark distributions at the same initial scale have been 
parametrized in line with Regge theory and the ideas 
leading to the gluon distribution. We have used DGLAP 
equations to evolve quark and gluon distributions to 
higher Q 2 values. The charm and proton structure func- 
tions computed in the small- X region are in good agree- 
ment with experimental HERA data over a large range 
of Q 2 values. 
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I. INTRODUCTION 

The understanding of deep inelastic scattering (DIS) in 
the small- a: regime remains one of the challenges in quan- 
tum chromodynamics (QCD). Perturbative and nonper- 
turbative physics are important for a complete picture of 
the small- a: limit of structure functions. In this work, we 
combine perturbative and nonperturbative approaches to 
describe charm and proton structure functions (or quark 
and gluon distributions) in the small- a; region. 

The basic idea is the description of high-energy scatter- 
ing in QCD by studying the eikonalized interaction of en- 
ergetic partons with external color fields. In ref. [1], this 
semiclassical method has been compared with the parton 
model to extract the leading order (LO) and next-to- 
leading order (NLO) gluon distribution. The main idea 
is simple and goes back to Muller [2]. One calculates 
the gluon production cross section for a scalar "photon" , 
which directly couples to the gluon field, in an external 
color field. The calculation of one-gluon production leads 
to the LO gluon distribution and the one of two-gluon 
production to the NLO gluon distribution. The LO result 
turns out to be a constant, xg^(x, Q 2 ) oc const., char- 
acterising the averaged local field strength in the proton. 
The NLO result shows a logarithmic increase at small 
x, xg^\x,Q 2 ) oc ln(l/a;), and is sensitive to the large 
distance structure of the proton. 
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In order to make numerical estimates for the gluon dis- 
tributions, we relate the LO and NLO gluon distributions 
to the scattering of a dipole in the adjoint representation 
of SU (3) on a proton. Then, we use the loop-loop corre- 
lation model (LLCM) [3] to compute this adjoint dipole - 
proton cross section at the center of mass (cm.) energy 
°f v^o w 20 GeV. The two gluons emerging from the 
scalar "photon" represent the adjoint dipole and a fun- 
damental quark-diquark dipole models the proton in the 
LLCM. The correlation between the two dipoles given 
by Wegner- Wilson loops is evaluated in the Gaussian 
approximation of gluon field strengths. The perturba- 
tive interactions are obtained from perturbation theory 
and the nonperturbative gluon field strength correlator 
is parametrized in line with simulations o in the lattice 
QCD. The calculation of scattering cross sections in the 
loop-loop correlation model has been quite sucessful at 
low energies y/s^ « 20 GeV [3-7]. 

For small virtuality Q 2 , the perturbative wave function 
of the two gluons emerging from the scalar "photon" is 
unrealistically extended at the endpoints of the longitudi- 
nal momentum fraction a = 0,1. This problem is similar 
to the quark-antiquark wave function of the transverse 
photon. The phenomenological solution we propose is to 
give the gluon a constituent mass ma of the order of the 
rho mass. This mass modifies the wave function at small 
transverse gluon momenta k± ensuring "confinement" for 
the two gluons, but does not affect the perturbative part 
of the two-gluon wave function at high k± . 

We calculate the ^-dependence of the gluon distribu- 
tion at a scale Qq = 1.8 GeV 2 which corresponds to 
the upper limit where the nonperturbative input is still 
credible. The quark distributions at the same scale are 



parametrized in line with Regge theory [8]. We evolve 
both distributions to higher values of Q 2 using DGLAP 
equations [9]. With quark and gluon distributions, we 
compute charm and proton structure functions at small 
x for different Q 2 values in good agreement with experi- 
mental data. 

The outline of the paper is as follows: In section II 
we review the formulas for the calculation of the LO and 
NLO gluon distribution following [1], introduce a gluon 
mass into the formalism and rewrite the final expressions 
for gluon distributions in terms of adjoint dipole-proton 
cross sections. We compute the adjoint dipole-proton 
cross section in the loop-loop correlation model in sec- 
tion III. In section IV we determine the x dependence of 
the quark distributions at the initial scale Qq according 
to Regge theory. The intial gluon and quark distributions 
are evolved to higher virtualities Q 2 using the DGLAP 
equation. In section V we present the results for the 
charm and proton structure function versus Bjorken-a; at 
different Q 2 values. Finally, in section VI, we summarize 
our results. 



II. 



THE SEMI-CLASSICAL GLUON 
DISTRIBUTION 



In this section we review the formulas for the compu- 
tation of the LO and NLO gluon distribution in the semi- 
classical approach and parton model following the work 
by Dosch, Hebecker, Metz and Pirner [1]. In addition, 
we introduce a gluon mass in the formalism to take into 
account the spatial localization of the two gluons due to 
confinement. We derive the modifications due to massive 
gluons and rewrite the final results for the gluon distri- 
butions in terms of adjoint dipole-proton cross sections. 



A. Gluon distribution at leading order 

In order to extract the gluon distribution, we consider 
a scalar field \ ( or scalar "photon") coupled directly to 
the gluon field through the interaction lagrangian 



X 



X tr(F M „FH 





(a) 



(b) 



FIG. 1: The process \ — Y g in an external color field: 
(a) one gluon exchange, (b) resummed gluon exchange. 



A A = A b t b the external color field, C A = N c the L-asimir 
operator in the adjoint representation, N c the number 
of colors and e± the polarisation vector of the outgoing 
gluon. 

Resumming the gluon exchange to all orders as shown 
in Fig. 1(b), the gluon initially created at the vertex \gg 
aquires a non-abelian eikonal factor on the way through 
the external color field 



U; 



(oo,a; + ) ( 6 -L 



P exp 



F 

■'.14 



dx + A_(x + , bj 



The path ordering along the way x + is denoted by P. 
With 

W b A (r A ) = U A (b ± -r A /2)U A \b ± +r A /2) - II , (1) 

where the eikonal factors U A (b± — r A /2) = U A _ OQ (b± — 
r A /2) and U A ^ (b±+r A /2) come from the scattering of two 
gluons (moving in opposite directions) with transverse 
distance r A off the proton, the semiclassical (sc) cross- 
section for gluon production at leading order becomes 
(cf. ref. [1]) 



A 2 



4g 2 C A 



J d 2 b ± \[d rA W A (r A )] rA=0 \ 



(2) 



In the following an average over the external color fields 
underlying the quantity Wj^ is implicitly understood. 

On the other hand, the leading-order parton model 
(pm) cross section for the partonic process \9 —> 9 shown 
in Fig. 2 with the proton described by the gluon distri- 



with Fpj, = F° v t a , the gluon field strength F" = 
[d^Al - d v A% - gf aie A b M Al] , the gluon field A a and the 
SU(N C ) group generators in the adjoint representation 
t a . In the high-energy limit, the scattering of a x-particle 
off an external color field A a as shown in Fig. 1(a), has 
the following amplitude 

T a (M = ^jdx + tx [t a (e* ± d ± )A A (x + ,b ± )] • 

Here b± denotes the impact parameter in transverse 
space, q the four momentum vector of the x - P ar ticle, 
x + = xq + X3 the light-cone variable, d± = d/db±, 



FIG. 2: The process \g — > g in the parton model, 
bution is given by 

a^ = ^xg<°Hx,Q 2 ). 



(3) 



Identifying the semiclassical cross-section (2) with the 
partonic one (3), one finds the leading-order gluon dis- 
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tribution 

xg^(x,Q 2 ) 



1 1 

2ir 2 a s 2CU 



JdH ± \[d rA W b A ± (r A )] rA=0 \ 



(4) 

Using the expression for the adjoint dipole (d) - proton 
(p) cross section in the semiclassical approach [1, 11] 



°%{u) = -\ / d 2 b ± trW^(ru) 
and the identity 

J d 2 b ± \[d rA W* (r A )] rA=Q | 2 = - \dl A J d 2 b ± tr (r A , 
one can rewrite eq. (4) in the useful form 



(5) 



xg {0) (x,Q 2 ) 



4jr 2 a s 2C A 



>r A =0 



r A =0 

(6) 



(7) 



This equation shows that the LO gluon distribution de- 
pends on the small distance behavior of the adjoint 
dipole-proton cross section. To get numerical results for 
the gluon distribution, we use in the next section a model 
to compute the adjoint dipole-proton cross section. 



B. Gluon distribution at next-to- leading order 



One can show [1] that the amplitudes in Fig. 4 produce 
the following NLO contribution to the cross-section 



(x,Q 2 



X 



3-2(2tt)« / a(l-a)I dk ' l I 



tfbi 



(8) 



/ 



N 2 + k\ +m 



2 yv bj 



- k±) 



where a and 1— a are the longitudinal momentum carried 
by the two gluons, k± and k' ± are the transverse momenta 
of one of the two gluons before and after the interaction 
with the external field, N 2 = a(l - a)Q 2 and W£ (k±) 
is the Fourier transform of W b A (r A ) given in eq. (1). 

In equation (8) we have introduced a gluon mass ma 
into the gluon propagator. This mass mimics the con- 
finement of the two gluons emerging from the x-particle 
and modifies their perturbative wave function in the non- 
perturbative, small k± region for small virtualities Q 2 or 
momentum fractions a = 0,1. Phenomenologically, the 
gluon mass has a similar effect as a constituent quark 
mass in the perturbative quark-antiquark wave function 
of the transverse photon [10]. We use for the gluon mass 
mo = 770 MeV, i.e., approximately half of the glueball 
mass. 

To compute the NLO contribution to the gluon dis- 
tribution in the parton model, we have to calculate the 
diagrams shown in Fig. 5. With z = Q 2 /{Q 2 + s) = x/y, 



For the semiclassical calculation of the gluon distribu- 
tion at next-to-leading order, the three diagrams shown 
in Fig. 3 are relevant in the high energy limit. Resum- 
ming the interaction with the external field to all or- 
ders, i. e., repeating the step leading from Fig. 1(a) to 
Fig. 1(b), we obtain the result shown in Fig. 4. In the 
first diagram the incoming x-partide splits into two glu- 
ons before interacting with the target. The subsequent 
scattering of the two gluons off the exteral color field 
is treated in the eikonal approximation. In the second 
diagram the two fast gluons are created through a X999 
vertex in the space-time region of the external color field. 




FIG. 5: Diagrams contributing to the parton model at NLO. 

where y is the fraction of the target momentum carried 
by the struck gluon, the parton model cross-section at 
NLO reads 






FIG. 3: Diagrams contributing to the x 
the high-energy limit. 



gg amplitude in 



with a 



(0) 



S(l-x)+aW(x,Q 2 )\ yg b (y), 

(9) 

(7rA^)/4 and the d = (4 + e) dimensional 
coupling Arf = A/_i~ £ / 2 . Since we are interested in the 
high energy limit, we extract only the leading term in 
ln(l/a;) from the expression in (9). We regularize the 
bare gluon distribution gb{x) in the MS scheme 





FIG. 4: Relevant contributions in the NLO semiclassical cal- 
culation. 



g b {x) = g(x,fi ) 



+ je — ln(47r) 



g(y,n 2 



with P gg denoting the gluon-gluon splitting function. 
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The partem model cross-section at NLO finally becomes 

(10) 



inhere 



7 d24 x [ ^r9(y, 

Jx z 



6(1 -z) + 



2tt 



Q 2 



where 



CfHz) = P aa {z)\n 



l-z\ IIC A 



z J Qz(l — z) 



Expanding the full gluon distribution up to the order 
ln(l/x) 

xg(x,n 2 ) = xg w (x,fi 2 ) + xg (1) (x, /i 2 ), 
and identifying (8) with (10), we obtain 



a. r 1 



P gg{z) ln[^)+Cr(z) 



(ii) 



where g^ (x, fi 2 ) is given by (4). When evaluating equa- 
tion (11), we keep only \n(l/x) terms, i.e., we use 



C™(z) 



2Ca 

: 

Z 

2C A 



In 



1\ 11 

" 12 



Using the variable z = Q 2 /(Q 2 + M 2 ), with M 2 = 
k l2 /{a(l — a)), and introducing a scale k 2 such that 
^qcd <S k 2 <C Q 2 , the first contribution in (11) can 
be split into a hard and a soft component [1, 11] 



L)? I a(l-a)I dk ' l I 

I 



d 2 b A 



^ N 2 8ij + 2h:k 



^W A (k' ± -k A 



N 2 + k 2 ± + m 2 G 



2 f 1 dz 



+-[ - r dk'if(k'D 

K Jx z Jo 



(12) 



tr [w^(ru)W^(0] e ik '^- r ^H(r A y A 



H 



r 2 r r2 



[a 2 K (a)] \b 2 K (b) 



2{r A -r' A f 



- 1 



{[aK^a)] [bK^b)] 
+ [aK 1 (a)] [b 2 K (b)] 
+ [a 2 K (a)] [Slfi(S)]} 



with a = mar A and 6 = mar' A . Note that the gluon 
mass does not influence the hard part which is calculated 
at leading twist. The lower bound of the z integration, 
x < 2, is a kinematical limit ensuring that the invariant 
mass of the two produced gluons cannot be larger than 
the total center-of-mass energy available [1, 11]. 

Inserting (12) in (11), the ln 2 (l/a;) terms from the 
semiclassical and the parton calculations cancel, so that 

x^W) = dk'lf(k' 2 ) 



+ ^C A 
it 



xg 



(o) 



For k 2 = fi 2 exp(ll/12), the second term also drops out 
and one obtains 

xgW(x,f) = ±ln\±) J dk'lf(k' 2 ). (13) 

After the integration over k' 2 and r' A , the NLO correction 
to the gluon distribution becomes [1, 11] 

1 , fl\ f°° dr 2 2 
x jd 2 b ± tr[w b A ± (r A )w£(r A )] (14) 



xgW(x,ti 2 ) 



where 



»oV) 
F(z) 



4ew _27E 

K\{z) + zK Q {z)K,{z) + Z ^K 2 (z) 
Z ^K 2 2 {z) - zK,{z)K 2 {z) + K 2 (z) 



In the large- N c limit [12], one obtains the relation 
Jd 2 b ± tr [w b A ± (r A )w£(r A )] = -2 J d 2 b ± tr Wj* (r 



a 



which allows us to rewrite the NLO correction in terms 
of the dip ole- proton cross section (5) 



xg (1) (x, n 2 



In 



dr 



F(m G r A ) af c (r A ). 

(15) 



The NLO gluon distribution depends on the large dis- 
tance behavior of the adjoint dipole-proton cross section. 
To obtain numerical results for the LO and NLO gluon 
distribution, we calculate the dipole-proton cross section 
at some low cm. energy y/sr] within the loop-loop corre- 
lation model in the following section. 



III. DIPOLE-PROTON CROSS SECTION FROM 
THE LOOP-LOOP CORRELATION MODEL 

In this section we use the loop-loop correlation model 
(LLCM) [3] to compute the cross section of an adjoint 
dipole scattering off a proton, a^(r A , so), at the cm. 
energy ^/s^ = 20 GeV. Then, we insert the resulting 
dipole-proton cross section in eqs. (7) and (15) to calcu- 
late the LO and NLO gluon distributions. 

The loop-loop correlation model is based on the func- 
tional integral approach to high-energy scattering in 
the eikonal approximation [4-7]. Its central elements 
are gauge-invariant Wegner- Wilson loops. With a phe- 
nomenological energy dependence, we have shown that 
the LLCM allows a unified description of high-energy 
hadron-hadron, photon-hadron, and photon-photon re- 
actions [3, 13] and of static properties of hadrons [14] in 
agreement with experimental and lattice data. In this 
work, we do not need a phenomenological parametriza- 
tion for the energy dependence since this is generated by 
the formalism outlined in the sections before. 

In the framework of LLCM the adjoint dipole - proton 
cross section reads [3, 15] 

^lcm^.So) = 2 j d2b ^j ^7 J d 2 r:F dz q \i! p (z q ,r r )\ 2 

x (l-S AT (s ,b ± ,r A ,z q ,r :F )) . (16) 

Here the correlation of two light-like Wegner- Wilson 
loops 



S Ajr (so,bj_,f A ,z q ,f :F ) 



(w A [c g ]w^[c v ]) c 



w A [c g ]) G {wr[c q ] /c 



describes the elastic scattering of a dipole in the funda- 
mental (J-) with a dipole in the adjoint (^4) representa- 
tion of SU(N C ). In the present case where a x-particle 
scatters off a proton, the adjoint color-dipole is given by 
the two gluons in the color-singlet state emerging from 
the x-partide and the fundamental color-dipole is given 
in a simplified picture by a quark and diquark in the 



proton. Each color-dipole is represented by a light-like 
Wegner- Wilson loop [16] 



W[C\ = Tr r Vex-p 



-ig£dz»Al(z)t a r 



where the subscript r indicates the representation of 
SU(N C ), Tr r = Tr r (---)/Trll I , is the normalized trace 
in the corresponding color-space with unit element 
and A lt (z) = A^(z)t^: represents the gluon field with the 
SU (N c ) group generators in the corresponding represen- 
tation, tf, that demand the path ordering indicated by 
V on the closed path C in space-time. Physically, the 
Wegner- Wilson loop represents the phase factor acquired 
by a color-charge in the SU (N c ) representation r along 
the light-like trajectory C in the gluon background field. 

The color-dipoles have transverse size and orientation 
r AT . The longitudinal momentum fraction of the dipole 
carried by the quark (gluon) is z q (z g ). The impact pa- 
rameter between the dipoles is [17] 

b± — Tg + (1 — Zg)r A — r q — (1 — z 9 )t> = r Acm — *> cra , 

where f q , f qq , f g , and r g are the transverse positions of 
the quark, diquark, gluon and the gluon moving in the 
opposite direction, respectively. With r> = r qq — r q and 
r A = r g — f g , the center of light-cone momenta of the 
two dipoles are given by r Acm = z g r g + (1 — z g )r g and 

cm = z q r q + (1 — z q )r qq . Figure 6 illustrates the (a) 
space-time and (b) transverse arrangement of the dipoles. 

The QCD vacuum expectation value (. . .)g in the loop- 
loop correlation function (17) represents functional inte- 
grals [5] over gluon field configurations: the functional 
integration over the fermion fields has already been car- 
ried out as indicated by the subscript G. The model we 
use for the QCD vacuum describes only gluon dynam- 
ics and, thus, implies the quenched approximation that 
does not allow string breaking through dynamical quark- 
antiquark production. 

The r> and z q distribution of the fundamental color- 
dipole in the proton is given by the proton wave function 
i[) p . We use for the proton wave function the phenomeno- 
logical Gaussian Wirbel-Stech-Bauer ansatz [18] 



* 9 (1 



^S 2 N P 



The constant N p is fixed by the normalization of the wave 
function to unity, the extension parameter S p is approxi- 
mately given by the electromagnetic radius of the proton 
and the width Az p = w/(^/2m p ) [18] is determined by 
(17) the proton mass m p and the value w = 0.35 — 0.5 GeV 



extracted from experimental data. We adopt the values 
Az p = 0.3 and S p = 0.86 fm which have allowed a good 
description of many high-energy scattering data in our 
previous work [3]. 

The computation of the loop-loop correlation for one 
dipole in the fundamental and the other one in the adjoint 
representation of SU(Nc) can be found in detail in [15]. 
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FIG. 6: High-energy scattering of a fundamental dipole with 
an adjoint dipole in the eikonal approximation represented 
by Wegner- Wilson loops: (a) space-time and (b) transverse 
arrangement of the Wegner- Wilson loops. The shaded areas 
represent the strings extending from the quark (gluon) to the 
diquark (antigluon) path in each color dipole. The thin tube 
allows to compare the field strengths in surface S q with the 
field strengths in surface S g . The impact parameter b± con- 
nects the centers of light-cone momenta of the dipoles. 



The main steeps of this computation are the transfor- 
mation of the line integrals in the Wegner- Wilson loops 
into surface integrals with the non-Abelian Stokes the- 
orem, a matrix cumulant expansion, and the Gaussian 
approximation of the functional integrals in the gluon 
field strengths. These steeps lead to the following result 



1 



{w A [c B ]) {w*[c t 



exp 



+ 
+ 



r. Jv c 
p[-4 



N c + 2 

2(N c +r 



which, for N c = 3, corresponds to the well-known SU(3) 
decomposition 

3®8 = 3ffil5 06. 

The function % := X F + Xnc + X™ F nas the following 
form [3]: 

X P =\g 2 (r g -r q )iD'^ (f g -f q ) 



^i=t z ^n'f 2 ) 1,7? 



9 ifg - fgg) iD^ (f g ~ Tgg) 

9 2 (r- g - fg) iD'^ ] (f- g - r 9 )] , 



(19) 



N p _ 7T 2 G 2 (1-K) 

inc ~~ 3(iV| - 1) 



- iD [ (2 \fg-r qq )-iD[ {2 \fg-f q )\ : 

(20) 



,NP 



(r A ■ »>) 



3(7V 2 - 1) 

x / dv A I dv T iD^ (r g 
Jo Jo 



(21) 



+ v A r A -r q - v T r T ) 



with the perturbative (iDp ') and nonperturbative 

(iD 1 ^ and iD^(z±)) correlation functions in transverse 
space 



D'^\z ± ) = ^-K (m G \z ± \), 



1 



2tt 



(22) 



^ (2) (4) =™ 4 [3 + 3(|4|/a) + (|4|/a) 2 ] exp(-|4|/a) 



iD^(z ± ) = 2na 2 [l + (\z ± \/a)]exp(-\z ± \/a) . 



(23) 



(24) 



We have introduced in the perturbative component \ p 
the same effective gluon mass ma = 0.77 GeV as before 
to limit the range of the perturbative interaction in the 
infrared region and a parameter M 2 = 1.04 GeV 2 which 
freezes the running coupling in the quenched approxima- 
tion at the value g 2 (z L )/(4n) = o> s (Jfc| = 0) = 0.4 [3], 



9 2 (z±) 



48tt 2 



(33 -2N } ) In [(\z ± \- 2 + M 2 )/A 2 QCD \ 



• (25) 



In Eq. (22) K denotes the th modified Bessel func- 
tion (McDonald function). The non-perturbative cor- 
relators (20) and (21) involve the gluon condensate 

G 2 := <^^(0)F^(0)> = 0.074 GeV 4 , the param- 
eter k = 0.74 that determines the relative weight of 
the two different components and the correlation length 
a = 0.302 fm that enters through the non-perturbative 
correlation functions D and D 1 . 
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The component \ F describes the perturbative interac- 
tion of the quark and diquark of the dipole in the proton 
with the two gluons of the adjoint dipole emerging from 
the x-particle as evident from (19) and Fig. 6b. The com- 
ponent Xnc nas the same structure as \ P ar >d gives the 
non-perturbative interaction between the quarks and glu- 
ons of the two dipoles. The component shows a dif- 
ferent structure: the integrations over v A and v T sum non- 
perturbative interactions between the strings (flux tubes) 
confining the quark and diquark or the two gluons in the 
dipoles as visualized in Fig. 6b. As shown in [14, 19, 20], 
the \c P component leads to color-confinement due to 
a flux tube formation between a static quark-antiquark 
pair. Manifestations of confinement in high-energy scat- 
tering have been analysed in [13]. 

Since the squared proton wave function \i[) p (z q , f>) | 2 
is invariant and the x _ f un ction changes sign under the 
replacement (r> — > — f>, z q — > 1 — z q ) 



x(b±,z g ,r A , 1 



-X(b±,z g , f A ,z q ,r% 



only the real part of the exponentials in eq. (18) survives 
in the integration over r T and z q so that one obtains 



W A [C g ]W^[C q 



1 



{w*[c B ]) {w*[c q \ 



+ 



N?-l 



N c + 2 
2{N C + 1) 

N c -2 
2(N C -1) 



2 



2 X 



(26) 



The above expression describes multiple gluonic interac- 
tions between two dipoles since (x P ) 2 represents the per- 
turbatively well-known two-gluon exchange and (x^) 2 
the non-perturbative two-point interaction in the dipole- 
dipole scattering [13]. The higher order terms in the ex- 
pansion of the cosine functions ensure the 5-matrix uni- 
tarity condition which becomes important at very high 
cm. energies [3, 15, 21]. 

The adjoint dipole-proton cross section obtained with 
the above ingredients at the cm. energy ^Jsq k, 20 GeV 
shows color-transparency for small dipole sizes, 



o d LLCM( s o,r A ) « 9.6 r 2 

and linear confining behavior at large dipole sizes, 

ct llcm( s o,^) oc \ f A \ . 

Inserting eqs. (27) and (25) into (7), one obtaines for the 
LO gluon distribution (7) at virtuality Ql = 1.8 GeV 2 

^ (0) ^^ = 4^^2^ 2 - 9 - 6 = °- 81 - < 29 > 

With our result for c*c Af ( s 0! , '-4)i the NLO gluon dis- 
tribution (15) at the same virtuality reads 



(27) 



(28) 



xg {1) (x,Q 2 ) = 0.89 In ( j- 



(30) 



IV. DGLAP EVOLUTION AT HIGH Q 2 

In this section we give the gluon and quark distribu- 
tions at an initial scale Ql Their evolution to higher 
values of Q 2 is obtained by the DGLAP equation 

h(x,Q 2 ) 




(31) 







with the splitting functions P xy being at leading order 
Q 2 independent. We use the resulting gluon and quark 
distributions to compute the charm and proton structure 
function at different x and Q 2 values. 

The gluon distribution computed in the previous sec- 
tion reads 

xg(x,Q 2 ) = A[l + B\n(l/x)] , (32) 

with A = 0.81 and B = 1.1 for Q 2 = 1.8 GeV 2 . Our re- 
sult, however, is only expected to be valid at low x values. 
For large x values, x > sgrv = 0.15, we use the Gluck- 
Reya-Vogt (GRV) gluon distribution [22]. To match to 
this gluon distribution at x = sgrv, we introduce a scale 
xq in our gluon distribution 

G(x, Ql) = xg(x, Q 2 ) = A [1 + B \n(x /x)] (33) 

which takes into account the neglected constant term in 
the NLO calculation where only the leading ln(l/a;) terms 
have been kept. For Qq = 1.8 GeV 2 and sqrv = 0.15, 
we obtain xq = 0.1454. 

To calculate the proton structure function F^(x,Q 2 ), 



Fi(x) 



e 2 [q(x) + q(x)] . 



flavours 



with the DGLAP evolution (31), we need quark and 
gluon distributions, since they are coupled to each other. 
For the computation of F$(x, Q 2 ), however, only two lin- 
ear combinations of quark distributions are required 

T = x(u + + c+ + f+) - x(d + +s + + b + ), 
E = x(u + +c+ + t + ) +x(d + +s + + &+), 

since 

5E + 3T 

F2 ' 

with q + = q + q and q = u,d, s, c, t, b. Performing linear 
combinations in (31), we can directly check that T, E and 
G evolve according to the following DGLAP equations 



Q 2 d Q 2T(x,Q 2 



2l f 1 



xd£ 



T(i,Q 



Q 2 d Q , 



e(*,q 2 ; 

G(x,Q 2 ] 



1 xdi (P l 



as f xdi 

aw* e 



2n s P q , 
p 



mQ 2 ) 

G(£,Q 2 ) 
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which means that T evolves alone, while E is coupled 
with G. 

We determine the s-dependence of the T and E dis- 
tribution at the scale Ql following the ideas from [23]. 
First, the T, E and G distributions are required to have 
a common singularity structure in the complex j-plane 
according to Regge theory. Since the initial gluon distri- 
bution (33) has a double pole in j = 1, consequently, also 
the T and E distribution are required to have a double- 
pole pomeron term. Secondly, we add a reggeon contribu- 
tion (coming from the exchange of meson trajectories a 
and /) to the quark distribution. In the gluon distribu- 
tion we negelect this term since the reggeon is expected to 
be constituted of quarks. Moreover, we expect that the 
pomeron, having vacuum quantum numbers, does not 
distinguish between quark flavours, i. e., the pomeron 
decouples from the T distribution. In contrast, the E 
distribution contains a pomeron and a reggeon compo- 
nent. The initial distributions, thus, read 



10 -4 10" 3 10" 2 



T{x,Ql) 
•Z(x,Ql) 



TX ao (l-x) a , 

[aln(l/x) + b + dx ao ](l-x) 



(34) X 



with the reggeon intercept cto = 0.4 and the powers a = 
3 and k = 2 of (1 — x) taking into account daughter 
trajectories in Regge theory. For x > sqrVj we again rely 
on the GRV distributions since our T and E distributions 
are only valid at small x. The parameters r and b are 
fixed to ensure continuity between our distributions and 
GRV's ones at x = »grv- Finally, we are left with 2 
parameters, a and d, which we determine by fitting the 
F%(x,Q 2 ) experimental data. 



V. RESULTS 

With the previous T, E and G distributions at Qg = 
1.8 GeV 2 as an initial condition in the DGLAP equation, 
we have fitted the F$ experimental data within the do- 
main 



f > Q 2 = 1.8GeV 2 
x < zqrv = 0.15 , 



(35) 



cos (6^ 



2xm T . 



> 



4 9 GeV 2 
2ml 



Here m p denotes the proton mass, 6 t is the scattering an- 
gle in the t-channel and cos(^) has been extended to the 
whole complex plane. The last condition of cos(^) taken 
from [24] ensures that Regge theory is applicable. We 
have used 922 experimental points coming from HI [25- 
27], ZEUS [28, 29], BCDMS [30], E665 [31] and NMC [32] 
to adjust our parameters within the given region. We 




FIG. 7: The charm structure function, _F|(x,Q 2 ), as a func- 
tion of the Bjorken- variable x at different virtualities Q 2 . 



have obtained the following values 

a = 0.45933 ±0.00425 , 

b = —1.9598 (fixed by continuity) , 

d = 6.0408 ± 0.0456 , 

t = 0.55628 (fixed by continuity) . 

In Figs. 8, 9 and 10 we show the results for the 
proton structure function, F^(x,Q 2 ), as a function of 
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FIG. 8: The proton structure function, F2(x,Q 2 ), as a function of the Bjorken-variable x at low virtualities Q 2 



the Bjorken-variable x between the virtualities Q 2 = 
1.8 GeV 2 and Q 2 = 3000 GeV 2 . Our results are in good 
agreement with experimental data for all x and Q 2 val- 
ues. 

To test our gluon distribution, we compute the charm 
structure function which depends on the gluon distribu- 
tion as follows [33] 



FC(x,Q 2 ) = 2e: 



f 



a s (Q 2 +4m 2 c ) 
2n 

dZg(£,Q 2 + 4m 2 c )f(x/^Q 2 



with 



[(4-/^(1 



x 



+ L [| - x 2 (l -x)+ fix 2 (l - 3x) - [i 2 z 3 \ 



2ml 

11 = -w 



2x/j 
\-x 
1 + v 



log , 

1 — V 

a = l + 2/i . 

We have adopted a value of 1.25 GeV for the charm 
quark mass. The predictions for F£ (x, Q 2 ) obtained from 



10 




10" 3 10" 2 10" 1 10" 3 10" 2 10" 1 10" 3 10" 2 10" 1 




X 



FIG. 10: The proton structure function, F% (x,Q 2 ), as a function of the Bjorken-variable x at high virtualities Q 2 . 
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our model are presented together with the experimental 
HERA data [28, 34, 35] in Fig. 7. The good agreement 
of our predictions for F^XjQ 2 ) with experimental data 
clearly shows that we obtain a resonable gluon distribu- 
tion. 



VI. CONCLUSIONS 

We have considered the interaction of a scalar "pho- 
ton", which directly couples to the gluons, with an exter- 
nal color field to extract the leading and next-to-leading 
order gluon distribution. We have closely followed the 
previous work by Dosch, Hebecker, Metz and Pirner [1], 
in which the semiclassical approach (where the partonic 
fluctuations of the "photon" interact with the proton in 
the eikonal approximation) has been compared with the 
parton model to get the gluon distribution of the proton. 
The leading order result is a constant while the next-to- 
leading order result has a ln(l/a;) rise at small x. 

We have been able to relate the leading and next-to- 
leading order gluon distribution to the short and long 
distance behavior of the cross section of a dipole in the 
adjoint representation of SU(3) scattering off a proton, 
respectively. In addition, a gluon mass has been intro- 
duced to take into account nonperturbative effects in the 



small-fci region of the perturbatively derived wave func- 
tion of the two gluons emerging from the scalar "photon" , 
in analogy to the constituent quark mass in the quark- 
antiquark wave function. We have computed the adjoint 
dipole-proton cross section in the loop-loop correlation 
model [3] to obtain numerical results for the gluon dis- 
tributions at the initial scale Qf, = 1.8 GeV 2 . 

Quark distributions at the same initial scale have been 
parametrized in line with Regge theory and the ideas 
leading to the gluon distribution. We have used DGLAP 
equations to evolve quark and gluon distributions to 
higher Q 2 values. The charm and proton structure func- 
tions computed in the small-s region are in good agree- 
ment with experimental HERA data over a large range 
of Q 2 values. 



Acknowledgments 

We would like to thank Arthur Hebecker and Stephane 
Munier for stimulating discussions. This research is par- 
tially funded by the INTAS project "Non-Perturbative 
QCD" and the European TMR Contract HPRN-CT- 
2000-00130. G. S. is supported by the National Fund 
for Scientific Research (FNRS), Belgium. 



[1] H. G. Dosch, A. Hebecker, A. Metz and H. J. Pirner, 

Nucl. Phys. B 568 (2000) 287. 
[2] A. H. Mueller, Nucl. Phys. B 335 (1990) 115. 
[3] A. I. Shoshi, F. D. Steffen and H. J. Pirner, Nucl. Phys. 

A 709 (2002) 131. 
[4] O. Nachtmann, Annals Phys. 209 (1991) 436. 
[5] O. Nachtmann, in "Perturbative and Nonperturbative 

Aspects of Quantum Field Theory," edited by H. Latal 

and W Schweiger (Springer Verlag, Berlin, Heidelberg 

1997) [hep-ph/9609365]. 
[6] A. Kramer and H. G. Dosch, Phys. Lett. B 252 (1990) 

669. 

[7] H. G. Dosch, E. Ferreira and A. Kramer, Phys. Rev. D 
50 (1994) 1992. 

[8] T. Rcgge, Nuovo Cim. 14, 951 (1959); 
T. Regge, Nuovo Cim. 18, 947 (I960); 
See for a modern overview of Regge theory and diffrac- 
tion the books Pomeron Physics and QCD by S. Don- 
nachie, G. Dosch, P. Landshoff and O. Nachtmann (Cam- 
bridge University Press, Cambridge, 2002), and High- 
Energy Particle Diffraction by V. Barone and E. Predazzi 
(Springer, Berlin, Heidelberg, 2002). 

[9] V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15 (1972) 781; 
L. N. Lipatov, Sov. J. Nucl. Phys. 20 (1975) 94; 
G. Altarelli and G. Parisi, Nucl. Phys. B 126 (1977) 298; 
Y. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 
[10] H. G. Dosch, T. Gousset and H. J. Pirner, Phys. Rev. D 

57 (1998) 1666. 
[11] M. Haas, "Berechnung der Gluonverteilung des Pro- 
tons unter Benutzung zweier nichtperturbativer Mod- 
elle," Diploma Thesis, Heidelberg 2000. 



[12] W. Buchmuller, T. Gehrmann and A. Hebecker, Nucl. 
Phys. B 537 (1999) 477. 

[13] A. I. Shoshi, F. D. Steffen, H. G. Dosch and H. J. Pirner, 
Phys. Rev. D 66 (2002) 094019. 

[14] A. I. Shoshi, F. D. Steffen, H. G. Dosch and H. J. Pirner, 
"Confining QCD strings, Casimir scaling, and a Eu- 
clidean approach to high-energy scattering," arXiv:hep- 
ph/0211287. 

[15] A. I. Shoshi, "Anatomy of QCD strings and saturation 

effects in high-energy scattering," arXiv:hep-ph/0302179. 
[16] F. J. Wegner, J. Math. Phys. 12 (1971) 2259; 

K. G. Wilson, Phys. Rev. D 10 (1974) 2445. 
[17] H. G. Dosch, T. Gousset, G. Kulzinger and H. J. Pirner, 

Phys. Rev. D 55 (1997) 2602. 
[18] M. Wirbel, B. Stech and M. Bauer, Z. Phys. C 29 (1985) 

637. 

[19] L. Del Debbio, A. Di Giacomo and Y. A. Simonov, Phys. 
Lett. B 332 (1994) 111. 

[20] M. Rueter and H. G. Dosch, Z. Phys. C 66 (1995) 245. 

[21] A. I. Shoshi and F. D. Steffen, "Saturation effects in 
hadronic cross sections," arXiv:hep-ph/0212070; 
A. I. Shoshi, F. D. Steffen and H. J. Pirner, "Gluon sat- 
uration and S-matrix unitarity," arXiv:hep-ph/0205343. 

[22] M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 
(1998). 

[23] G. Soyez, Phys. Rev. D 67 (2003) 076001. 

[24] J. R. Cudell, E. Martynov and G. Soyez, "t-channel uni- 
tarity and photon cross sections," arXiv:hep-ph/0207196. 

[25] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 13, 
609 (2000). 

[26] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 19, 



13 



269 (2001). 

[27] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 21, 
33 (2001). 

[28] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 
12, 35 (2000). 

[29] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. 

C 21, 443 (2001). 
[30] A. C. Benvenuti et al. [BCDMS Collaboration], Phys. 

Lett. B 223, 485 (1989). 
[31] M. R. Adams et al. [E665 Collaboration], Phys. Rev. D 



54, 3006 (1996). 
[32] NMC Collaboration: M. Arneodo et al, Nucl. Phys. 

B483 (1997) 3; Nucl. Phys. B487 (1997) 3. 
[33] R.G. Roberts, The Structure of the Proton, Cambridge 

University Press (1990). 
[34] C. Adloff et al. [HI Collaboration], Z. Phys. C 72, 593 

(1996). 

[35] J. Breitweg et al. [ZEUS Collaboration], Phys. Lett. B 
407, 402 (1997). 



